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The  solution  to  the  problem  of  the  interaction  of  a  spherical  wave  with  a  homogeneous  and 
isotropic  solid  plate  of  infinite  lateral  extent,  but  finite  thickness,  is  considered  theoretically. 
Both  the  source  and  the  plate  are  immersed  in  an  infinite,  inviscid  fluid.  Appropriate  boundary 
conditions  are  imposed  on  the  full  three-dimensional  elasticity  equations.  The  solution  is 
evaluated  numerically  for  a  variety  of  materials  for  a  1-kHz  incident  spherical  wave  and  for  a 
5-kHz  incident  spherical  wave.  For  the  numerical  study,  the  fluid  medium  is  taken  to  be  water. 
Under  certain  conditions,  “overpressures”  are  predicted  for  both  the  reflected  and  transmitted 
fields  (i.e.,  the  amplitude  of  the  reflected  pressure  and/or  the  transmitted  pressure  can  exceed 
the  maximum  value  of  the  amplitude  of  the  incident  pressure  on  the  plate  surface).  These 
overpressures  are  consistent  with  the  law  of  conservation  of  energy  in  the  sense  that,  for  a  plate 
composed  of  lossless  material,  the  total  incident  power  is  found  to  be  equal  to  the  sum  of  the 
total  reflected  power  plus  the  total  transmitted  power.  An  important  conclusion  of  this 
research  is  that  the  practice  of  attempting  to  reduce  the  influence  of  edge  diffraction  in  panel 
tests  by  using  samples  of  increasingly  larger  lateral  extent  may  result  in  measurements  that  are 
substantially  corrupted  by  wave-front  curvature  effects,  particularly  if  the  sample  panel 


includes  a  steel  backing  plate. 

PACS  numbers:  43.20.Fn,  43.30.Gv 

LIST  OF  SYMBOLS 

Pr*Pt 

reflected  pressure  and  transmitted  pres¬ 

A 

vector  potential  in  the  solid 

P' 

sure,  respectively 

a  general  pressure  variable  equal  to 

^0 

jPae-^z' 

Pi  +  Pr  when  evaluated  at  z  =  0,  and 

a(0)MP),c(P), 

d(P),e(p)M) 

expansion  functions 

Po 

equal  to  p,  when  evaluated  at  z  =  —  /. 
incident  spherical  wave  pressure  ampli¬ 

c„c , 

longitudinal  speed  and  shear  speed  in  the 

tude  coefficient;  p,  =  P<£'e~ik*e>k'/r 

solid,  respectively 

r,&rt 

cylindrical  polar  coordinates  (see  Fig.  1 ) 

E„Er,E, 

incident  power,  reflected  power,  and 

/ 

radius  of  the  “circle  of  observation”  (see 

j 

transmitted  power,  respectively 
imaginary  unit  =  V  —  1 

R 

text  for  definition) 

reflection  coefficient  evaluated  at  the  ori¬ 

J0(z)J^z) 

cylindrical  Bessel  functions  of  the  first 

gin  of  coordinates 

kind  of  argument  z  and  orders  0  and  1, 

T 

transmission  coefficient  evaluated  at  the 

k 

respectively 

real  wavenumber  in  the  inviscid  fluid  sur¬ 

u 

point  (0,0,  —  /) 

particle  displacement  in  the  solid 

rounding  the  plate 

z' 

z  coordinate  of  the  spherical-wave  source 

k„k, 

complex  longitudinal  wavenumber  and 

Aji 

Lame  parameters 

complex  shear  wavenumber  in  the  solid, 

Pi 

density  of  the  solid 

^2 

respectively 

wavenumber  in  the  fluid  slab 

&r>00>&z* 
~ez>ri.  >TrO 

stress  components 

l 

plate  thickness 

+ 

scalar  potential  in  the  solid 

Pi 

incident  pressure 

0) 

angular  frequency  of  the  spherical  wave 

INTRODUCTION 

Generally,  the  plane-wave  reflection  coefficient  and  the 
plane-wave  transmission  coetficient  are  the  quantities  of  in¬ 
terest  in  panel  measurements.  Plane-wave  conditions  are  of¬ 


ten  established  by  utilizing  the  farfield1  of  a  source.  How¬ 
ever,  due  to  the  limited  size  of  test  facilities,  it  is  not  always 
practical  to  achieve  the  desired  separation  between  the 
source  and  the  test  panel.  Therefore,  it  is  important  to  deter- 
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mine  the  potential  ramifications  of  performing  panel  mea¬ 
surements  in  the  nearfield'  of  a  source.  One  of  the  simplest 
cases  to  consider  is  that  of  the  sound  radiation  from  a  point 
source  interacting  with  a  single-layer,  homogeneous,  iso¬ 
tropic  solid  plate  of  infinite  lateral  extent  but  finite  thick¬ 
ness. 

The  interaction  of  a  spherical  wave  with  a  plane  has  long 
been  of  interest  to  the  geophysics  community.2-6  However, 
the  results  of  this  previous  work  are  not  immediately  appli¬ 
cable  to  panel  measurements  for  a  variety  of  reasons.  Often, 
only  high-frequency  or  farfield  limiting  cases  have  been  con¬ 
sidered.  Much  of  this  previous  work  has  not  included  the 
influence  of  shear  waves  in  the  solid,  or  has  neglected  loss,  or 
has  considered  only  semi-infinite  solid  media.  Rarely  has  a 
transmitted  wave  been  computed,  and  generally  not  for  the 
case  in  which  the  medium  behind  the  solid  (with  respect  to 
the  source)  is  also  an  inviscid  fluid.  In  order  to  be  useful  in 
evaluating  panel  measurements,  one  needs  the  general  solu¬ 
tion  to  the  problem  of  a  spherical  wave  interacting  with  a 
lossy,  finite-thickness  solid  that  is  both  homogeneous  and 
isotropic.  The  source  and  the  solid  must  both  be  immersed  in 
an  inviscid  fluid  medium.  It  is  the  purpose  of  the  present 
article  to  discuss  the  consequences  of  the  solution  to  this 
problem  for  panel  measurements. 

The  approach  taken  to  obtain  the  solution  is  to  impose 
the  appropriate  boundary  conditions  on  the  full  three-di¬ 
mensional  elasticity  equations.  Solutions  to  elasticity  prob¬ 
lems  are  rarely  based  on  the  complete  elasticity  theory, 
owing  to  the  complexity  of  the  attendant  mathematics.  How¬ 
ever,  due  to  the  relatively  recent  advent  of  computer  pro¬ 
grams  that  are  capable  of  directly  manipulating  mathemat¬ 
ical  symbols,  this  difficulty  is  no  longer  insuperable.  Unlike 
computer  languages  like  fortran  and  basic,  such  symbol 
manipulation  programs  do  not  require  numerical  values  to 
be  assigned  to  the  symbolic  quantities  of  interest.  Thus  one 
can  bring  the  full  power,  speed,  and  memory  capabilities  of  a 
modem  computer  to  bear  on  problems  involving  the  manipu¬ 
lation  of  quantities  containing  thousands  of  algebraic  sym¬ 
bols.  In  addition  to  being  capable  of  performing  the  standard 
algebraic  operations  of  expansion,  factorization,  rationali¬ 
zation,  etc.,  many  such  programs  can  also  perform  several 
calculus-type  operations,  such  as  differentiation  (of  very 
complicated  expressions)  and  integration  (of  somewhat  ele¬ 
mentary  expressions).  Of  course,  such  programs  are  really 
only  “slaves”  that  can  perform  calculations  that  are  capable 
of  algorithmic  description,  such  as  computing  the  derivative 
of  an  expression  containing,  say,  100  factors  (a  task  that  may 
be  beyond  human  capability).  However,  such  programs  are 
incapable  of  performing  tasks  that  require  human  talent, 
such  as  choosing  an  appropriate  series  expansion  to  satisfy 
both  a  given  differential  equation  and  the  appropriate 
boundary  vCuditions.  They  are,  of  course,  completely  in¬ 
capable  of  choosing  which  differential  equation  to  solve,  as 
weii  as  which  boundary  conditions  are  appropriate.  Thus 
such  symbol  manipulators  bear  approximately  the  same  re¬ 
lationship  to  the  process  of  obtaining  the  solution  to  a  phys¬ 
ical  problem  as  word  processing  programs  bear  to  the  pro¬ 
cess  of  creating  a  manuscript. 

The  theoretical  results  discussed  in  this  article  were  de¬ 


rived  with  the  help  of  the  SMP™  computer  program.  This 
program  also  has  the  added  benefit  of  being  able  to  directly 
generate  fortran  subroutines  that  can  numerically  evalu¬ 
ate  the  resulting  complicated  symbolic  expressions.  Thus, 
once  the  symbolic  expressions  for  the  boundary  conditions 
and  the  assumed  solutions  have  been  defined,  the  process  of 
imposing  the  boundary  conditions  on  the  assumed  solutions 
is  relatively  straightforward.  However,  the  problem  of  nu¬ 
merically  evaluating  the  required  integrals  must  still  be 
tackled  by  more  conventional  means. 

The  theory  is  developed  in  Sec.  I.  Numerical  results  are 
presented  in  Sec.  II.  Some  consequences  of  the  theoretical 
results  for  panel  measurements  are  given  in  Sec.  Ill  Finally, 
a  summary  and  the  conclusions  are  given  in  Sec.  IV. 

I.  THEORY 

The  theoretical  approach  used  here  is  the  same  as  that 
used  in  Ref.  1,  in  which  the  scattering  of  a  spherical  wave  by 
an  elastic  cylinder  is  considered.  However,  unlike  that  case, 
the  differential  equations  in  the  present  problem  are  solved 
exactly.  Also,  unlike  Ref.  1,  the  numerical  investigation  con¬ 
sidered  here  includes  the  effects  of  nonzero  loss  in  both  the 
bulk  and  shear  moduli  of  the  material  composing  the  plate 

A.  Differential  equations 

The  geometry  of  the  scattering  problem  and  the  coordi¬ 
nate  system  are  presented  in  Fig.  1 .  Following  the  procedure 
of  Ref.  1,  we  begin  the  present  investigation  by  first  consider¬ 
ing  the  equation  of  motion  of  waves  in  a  solid  elastic  medi¬ 
um7 

(A  +  2*u)V(V-u)  -yuVx(Vxu)  (1) 

ot 

We  assume  that  the  displacement  u  in  the  solid  can  be  re' 


♦z 


Fi<u.  i.  A  ^uart  li  Tnd  *•  joiiu 

thickness  /,  with  infinite  lateral  dimensions.  An  infinite,  inviscid  fluid  medi¬ 
um  is  located  on  both  sides  of  the  plate.  A  spherical-wave  source  is  located 
at  an  offset  distance  z'  from  one  side  of  the  plate.  The  point  of  contact  of  the 
perpendicular  between  the  spherical-wave  source  and  the  plate  defines  the 
origin  of  coordinates.  The  coordinate  system  is  cylindrical  polar,  with  coor¬ 
dinate  variables  r.  0,  z.  For  clarity  in  defining  the  polar  angle,  a  rectangular 
coordinate  system  is  also  depicted. 
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solved  via  the  usual  Helmholtz  theorem  into  scalar  and  vec¬ 
tor  potentials  tf>  and  A,  respectively,  such  that 

u  =  —Vt/i  +  Vx  A.  (2) 


Here,  if>  and  A  are  required  to  be  solutions  of  the  differential 
equations 


vV- 


1  a24> 

c?  dt 2 


=  o 


(3) 


and 


v2a  1^-0. 
ej  dt2 

(4) 

We  can  rewrite  Eq.  (4)  by  expressing  V2A  in  component 
form,8  giving 

V2f  Ae  2  dA,  _ 

*  r2+r2dd  ' 

(5) 

Ar  2  dAg  _ 

WAr  r2  r2  d0 

(6) 

and 

V2A,  =  —k\A,. 

(7) 

Thus  the  coupled  set  of  four  partial  differential  equations  to 
be  solved  consists  of  Eqs.  (3)  and  (5)-(7).  Note  that  Eq.  (3) 
yields  the  Helmholtz  equation  if  harmonic  time  dependence 
is  assumed.  Note,  in  addition,  that  Eq.  ( 7 )  is  also  a  Helmholtz 
equation.  However,  Eqs.  (5)  and  (6)  arenor  Helmholtz  equa¬ 
tions.  One  consequence  of  this  fact  was  that  the  problem  of 
the  scattering  of  a  spherical  wave  by  an  elastic  cylinder, 
considered  in  Ref.  1,  was  only  solved  approximately. 

From  the  symmetry  of  Fig.  1,  it  is  clear  that  none  of  the 
solutions  to  the  present  problem  will  depend  on  the  coordi¬ 
nate  variable  0.  Hence,  the  terms  containing  a  derivative 
with  respect  to  this  variable  in  Eqs.  ( 5 )  and  ( 6 )  are  both  zero. 
This  means  that  these  equations  simplify  to 

V2Ag-Ag/r=-k2sAg  (8) 

and 


V2A ,  —  Ar/r  =  -k]Ar 


(9) 


respectively.  Thus,  unlike  the  analogous  situation  arising  in 
the  case  of  the  cylindrical  scatter,  the  differential  equations 
for  the  functions  A0  and  A,  in  the  present  problem  are  not 
coupled.  However,  Eqs.  (8)  and  (9)  still  are  not  Helmholtz 
equations,  due  to  the  presence  of  the  —  A  /r2  terms  in  each 
equation.  We  will  therefore  follow  the  usual  procedure  used 
to  solve  such  equations;  namely,  we  will  introduce  the  as¬ 
sumption  of  separation-of-variables.  Since  Eqs.  (8)  and  (9) 
are  of  identical  form,  we  need  to  consider  in  detail  only  one 
of  these  equations.  InEq.  ( 8 ) ,  we  assume  that  the  function  /l 
is  separable  into  functions  R(r),  9(0),  and  Z(r)  such  that 
Ag  =  R(r)Q(9)Z(z).  Substitution  of  this  assumed  form 
into  Eq.  (8)  produces  the  separated  set  of  equations 


1  d2Z(z) 
Z(r)  dz 2 


(10) 

_ I 


Pr  = 


<*(0)Jo(0r)exp[  -jyjk’1  —fil(z'  -  z)]0d/3 


1  rf26(0)  _ 

0(0)  d9 2 


(11) 


1  d(dR{r)\ 
rR(r)dA  dr  ) 


-l=-*s2,  (12) 


where  m 2  and  k  2  are  <*paration  constants.  Equations  (10) 
and  (11)  admit  complex  exponential  and  trigonometric  so¬ 
lutions,  where  the  separation  constant  m  is  required  to  be  an 
integer  by  the  periodicity  condition  associated  with  the  coor¬ 
dinate  variable  0.  Equation  (12)  admits  solutions  that  are 
Bessel  functions  of  irrational  order  v'  l  +  m1.  However,  due 
to  the  assumed  independence  of  the  solutions  from  the  vari¬ 
able  0,  the  separation  constant  m  must  vanish.  If  we  intro¬ 
duce  a  new  constant  k 2  such  that 

k2w=k2s-k\,  (13) 


the  radial  equation  simplifies  to 


'-Mr 

r  dr\ 


*7r)+{k‘--7hr>-°- 


(14) 


The  only  solution  to  Eq.  ( 14)  that  obeys  the  physical  require¬ 
ment  that  R (r)  remain  finite  for  r  =  0  (see  Fig.  1 )  is 

R(r)  =CJ,(ku,r),  (15) 


where  C  is  an  arbitrary  constant.  This  means  that  the  solu¬ 
tions  for  the  functions  Ag  and  Ar  can  be  built  up  out  of  com¬ 
plex  exponential  (or  trigonometric)  functions  and  the  func¬ 
tion  This  matter  will  be  considered  further  in  the  next 
subsection. 


B.  Assumed  forms  for  the  solutions 

In  order  to  be  able  to  solve  the  differential  equations 
represented  by  Eq.  (3)  ( under  the  assumption  of  harmonic 
time  dependence)  and  Eqs.  (7)-(9),  it  is  first  worthwhile  to 
choose  an  appropriate  expansion  for  the  incident  spherical 
wave.  This  is9 


/0 (/?r) exp  [/V&  2  -  /? 2  (z'  -  z)  ] 

4k^W 


P  dP, 

(16) 


which  is  valid  for  z<z‘. 

In  constructing  the  solutions  to  the  present  problem,  it 
is  important  to  note  that  this  process  is  greatly  simplified  in 
the  present  case  since  all  such  solutions  can  be  expressed  in 
terms  of  standard  wavefunctions.  It  is  also  important  to  take 
cognizance  of  the  symmetry  of  the  problem.  In  the  present 
case,  symmetry  requires  the  solutions  to  be  independent  of 
the  0  coordinate.  Also,  it  is  important  to  bear  in  mind  the 
boundary  conditions  that  must  be  satisfied.  These  will  be 
considered  in  the  next  subsection.  We  will  therefore  proceed 
to  assume  that  the  solutions  to  the  present  problem  may  be 
expressed  in  the  form 


(17) 
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(17) 


P'=  J[ 


__  f"  a(^)J0(j9r)exp[  -yVfe7  -fF(z' -z)]0d0 
rp—p? 

_  f°°  y0(jgr){c(jtf)exp[y/A:^-^-(2'-2)]  +  rf(ff)exp[  f  -  /?2(z'  -  z)  ] 

A  _  p  ■/|(^r){g(j9)exp[yvTi2  -/g'(z'  -z)  ]  +/(ff)exp[  -/V*? -/92(z' -z)]}*3  4S 

Jo  ’ 


(19) 

(20) 


and  A2  —  Ar  =  0.  Here,  a,  b,  c,  d,  e,  and  / are  six  unknown 
functions  to  be  determined. 

Of  course,  in  order  to  justify  these  assumed  forms,  two 
conditions  must  be  satisfied,  hirst,  they  must  represent  par¬ 
ticular  solutions  of  the  partial  differential  equations.  Sec¬ 
ond,  they  must  be  capable  of  satisfying  all  the  boundary  con¬ 
ditions  involved  in  the  physical  problem.  We  take  up  the 
latter  point  in  the  next  subsection. 

To  see  that  the  given  assumed  forms  do  in  fact  satisfy 
the  differential  equations,  note  first  that  the  quantities  pr 
and  p,  represent  pressure  waves  in  the  fluid  surrounding  the 
solid.  Therefore,  the  differential  equation  that  governs  the 
behavior  of  these  functions  is  the  Helmholtz  equation. 
Equations  (17)  and  (18)  areclearly  particularsolutionsof  the 
Helmholtz  equation  because  they  are  constructed  from  stan¬ 
dard  wavefunctions.  [Note  also  the  similarity  of  Eqs.  (17) 
and  (18)  to  the  expansion  for/i,  given  inEq.  (16).]  The  same 
argument  applies  to  the  expansion  for  the  function  \p,  given 
in  Eq.  ( 1 9 ) ,  since  this  function’s  governing  differential  equa¬ 
tion,  Eq.  (3),  is  also  a  Helmholtz  equation,  under  the  as¬ 
sumed  harmonic  time  dependence.  The  solution  for  Ae  has 
been  constructed  to  be  consistent  with  the  discussion  given 
inconnection  with  Eqs.  ( 10)— ( 15).  Hence,  Eq.  (20)  is  seen  to 
be  a  particular  solution  of  Eq.  (8).  Finally,  note  that  the 
equations/*,  =Ar  =  Oare  the  trivial  solutions  to  Eq.  (7)and 
Eq.  (9),  respectively.  Thus  it  only  remains  to  be  shown  that 
these  assumed  forms  can  satisfy  the  boundary  conditions  of 
the  problem.  We  consider  this  matter  next. 


C.  Boundary  conditions 


There  are  three  boundary  conditions  for  the  present 
problem,  each  of  which  must  be  imposed  at  each  of  the  two 
fluid-solid  interfaces.  [This  explains  why  six  unknown  ex¬ 
pansion  functions  a-f  are  needed  to  formulate  the  assumed 
solutions  represented  by  Eqs.  ( 17)— (20). ]  These  three 
boundary  conditions  are:  (i)  the  continuity  of  normal  dis¬ 
placement,  (ii)  the  continuity  of  normal  stress,  and  (iii)  the 
vanishing  of  shear  stress.  By  combining  the  general  expres¬ 
sion  for  strain  written  in  cylindrical  coordinates10  with 
Hooke’s  law  for  isotropic  media,11  we  obtain  the  stress- 
strain  relationships 


dur 

o,  =  2/t  —  +  /iV*u, 
dr 

(21) 

lu 

+  rlV-U, 

(22) 

r  dd 

du 

<7,  =  2p  — -  +  /lV*u, 
dz 

(23) 
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To* 

(du2  dug\ 
Krd&  dz  ) 

(24) 

=/'( 

'du,  | 

.  dz  ot  J 

(25) 

and 

<  du„  ue 

1  du,\ 

(26) 

‘  r6  —  fM 

i.  dr  r 

r  dOJ 

It  is  next  helpful  to  reexpress  the  stress  components  in 
terms  of  the  potentials  A  and  tp.  This  may  be  done  with  the 
help  of  Eq.  (2).  Also,  it  should  be  noted  that  attempting  to 
directly  impose  the  boundary  conditions  will  result  in  a  set 
of  equations  that  is  difficult  to  solve.  A  more  convenient  set 
of  equations  is  obtained  if  the  equation  that  imposes  contin¬ 
uity  of  normal  displacement  is  differentiated  with  respect  to 
r.  ( This  is  a  legitimate  operation  because,  although  this  equa¬ 
tion  is  evaluated  at  a  particular  value  of  z,  it  must  hold  for  all 
values  of  r.  Hence,  each  side  of  the  equation  is  a  continuous¬ 
ly  differentiable  function  ofr.)  Equation  (8)  can  be  used  to 
replace  combinations  of  terms  in  the  boundary  condition 
equations  with  simplified  expressions.  If  we  also  utilize  the 
symmetry  requirement  that  the  solutions  be  independent  of 
0,  and  if  we  further  take  advantage  of  the  assumption  that 
Az  =Ar  =  0,  the  boundary  conditions  assume  the  form 


dzdr  V  dr  dz  dz~  / 

(for  the  continuity  of  normal  displacement), 


(for  the  continuity  of  normal  stress),  and 


0=/r( 


2  _£V_ 

dz  dr 


d2A. 

2 - - 

dz 2 


(27) 


(28) 


(29) 


(for  the  vanishing  of  shear  stress). 


Note  that  Eqs.  (27)-(29)  are  written  so  that  the  left-hand 
sides  refer  to  the  fluid  and  the  right-hand  sides  refer  to  the 
solid  [e.g.,  the  left-hand  side  ofEq.  (28)  represents  the  stress 
in  the  fluid,  while  the  right-hand  side  of  this  equation  repre¬ 
sents  the  normal  stress  in  the  solid  ] .  Each  of  these  boundary 
conditions  must  be  imposed  both  at  z  =  0  and  z  =  —  l  ( refer 
to  Fig.  1 ) .  The  integral  equations  that  result  from  substitut- 
ingEqs.  ( 16)-(20)intotheboundaryconditionsofEqs.  (27)- 
(29)  can  be  converted  into  algebraic  equations  by  using  the 
orthogonality  property  of  the  Bessel  function.  This  results  in 
a  set  of  six  algebraic  equations  in  the  six  unknown  functions 
a-f.  As  previously  remarked,  this  procedure  was  implemen¬ 
ted  with  the  help  of  the  symbolic  manipulation  computer 
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program  SMP.  This  program  was  also  used  to  analytically 
solve  the  resulting  algebraic  set,  and  to  generate  fortran 
subroutines  that  permit  numerical  evaluation  of  the  expan¬ 
sion  functions  a(0)  and  b(p). 

Note  that,  since  the  algebraic  set  of  six  equations  in  the 
six  functional  unknowns  a-f  is  capable  of  analytical  solu¬ 
tion,  this  completes  the  demonstration  that  the  assumed 
forms  of  the  solutions,  given  by  Eqs.  ( 17)-(20),areappropri- 
ate  for  the  present  problem.  The  resulting  analytical  solu¬ 
tions  for  the  expansion  functions  are  extremely  complicated 
and  unwieldy,  and  are  not  particularly  illuminating.  There¬ 
fore,  we  will  not  display  the  analytical  expressions,  and  will 
choose  instead  to  investigate  several  numerical  examples. 


II.  NUMERICAL  CALCULATIONS 

Preliminary  plots  of  the  real  and  imaginary  parts  of  the 
integrands  containing  the  expansion  functions  a{/3)  and 
b(0)  revealed  variations  [for  the  materials  and  frequency 
interval  of  interest  ( 1—5  kHz)  ]  that  are  similar  to  the  varia¬ 
tions  observed  in  the  expansion  functions  of  Ref.  1.  Namely, 
the  integrands  vary  slowly  over  most  of  the  integration  inter¬ 
val,  but  can  have  large  fluctuations  about  a  small  number  of 
points  (see  Fig.  3  of  Ref.  1).  Hence,  a  similar  numerical 
scheme  was  applicable  in  the  present  case  for  evaluating  the 
integrals  of  Eqs.  (17)  and  ( 18).  (These  equations  give  repre¬ 
sentations  of  the  reflected  pressure  and  the  transmitted  pres¬ 
sure,  respectively. )  In  summary,  the  numerical  scheme  in¬ 
volves  first  splitting  the  infinite  integration  interval  into  a 
finite  subinterval  plus  an  infinite  subinterval.  The  finite  su¬ 
binterval  is  then  searched  for  the  regions  of  most-rapid  vari¬ 
ation  of  the  integrand.  The  integral  over  the  finite  subinter¬ 
val  is  then  evaluated  using  a  Gauss-Legendre12  numerical 

_ I 


integration  procedure,  with  the  greatest  density  of  quadra¬ 
ture  points  being  located  in  the  regions  of  most-rapid  vari¬ 
ation  of  the  integrand.  The  small  contribution  to  the  total 
integral  from  *he  infinite  subinterval  is  then  evaluated  using 
a  15-point  Gauss-I  aguerre 13  numerical  integration  proce¬ 
dure.  See  Ref.  1  for  a  more  detailed  description  of  the  nu¬ 
merical  integration  process. 

As  was  the  case  for  the  solution  presented  for  the  expan¬ 
sion  function  in  Ref.  1,  the  solutions  for  the  expansion  func¬ 
tions  in  the  present  problem  are  extremely  complicated. 
However,  since  the  analytical  solutions  and  the  required 
fortran  subroutines  for  evaluating  them  numerically  were 
obtained  using  a  symbolic  computation  program,  the  likeli¬ 
hood  of  mathematical  slips  occurring  in  obtaining  the  solu¬ 
tions,  and  in  generating  the  required  computer  programs  to 
numerically  evaluate  the  results,  is  considerably  reduced. 
Nonetheless,  it  is  important  to  investigate  limiting  test  cases 
to  gain  confidence  in  the  reliability  of  the  software  before 
venturing  into  calculations  of  a  more  general  (and  more  in¬ 
teresting)  nature. 

A  very  useful  test  case  is  the  problem  of  the  interaction 
of  a  spherical  wave  with  a  finite-thickness  fluid  plane.  As  is 
the  case  for  the  solid,  an  infinite,  inviscid  fluid  medium  is 
located  on  both  sides  of  the  fluid  plane.  (The  two  fluids  in¬ 
volved  in  this  reduced  problem  are  of  course  assumed  to  be 
immiscible.)  The  solution  to  this  problem  is  relatively 
straightforward  to  obtain  by  analogy  to  the  solution  to  the 
related  problem  involving  an  incident  plane  wave.14 

The  functional  expansions  introduced  for  the  incident, 
reflected,  and  transmitted  waves  for  the  purpose  of  obtain¬ 
ing  the  solution  to  the  general  problem  [i.e.,Eqs.  ( 16)— ( 18), 
respectively]  are  also  applicable  to  the  finite-thickness  fluid 
plane  problem.  For  the  waves  internal  to  the  fluid  plane,  the 
following  expansion  is  also  useful15: 


P\nx 


Jo(.0r){g(P)e\p[jyjkl  -02(z' -z)\  +6(ff)exp[  —j^k\  -0~  (z'  -  z)  ]  }0  d0 
_  “  - 


(30) 


In  Eq.  ( 30) ,  k2  denotes  the  wavenumber  in  the  fluid  compos¬ 
ing  the  plane,  and  g  and  h  are  functions  to  be  determined.  [  Of 
course,  the  functions  a  and  b  of  Eqs.  (17)  and  (18)  are  also 
unknowns  that  must  be  redetermined  for  the  present  case.  ] 
The  boundary  conditions  for  the  finite-thickness  fluid  plane 
problem  are:  (i)  the  continuity  of  normal  displacement  and 
( ii )  the  continuity  of  pressure.  Each  of  these  conditions  must 
be  evaluated  both  at  z  =  0  and  at  z  =  —  /.  This  explains  why 
four  unknown  expansion  functions  a,b,g,  and  h  are  required 
in  the  representations  of  the  solutions  for  the  reduced  prob¬ 
lem. 

The  solution  to  the  fluid  plane  problem  described  above 
was  also  implemented  using  the  SMP  symbolic  computation 
program.  The  same  subroutine  that  was  used  to  numerically 
evaluate  the  integrals  required  to  compute  the  reflected- 
wave  pressure  and  the  transmitted-wave  pressure  in  the 
more  general  case  was  also  used  to  evaluate  the  associated 
integrals  arising  in  the  reduced  problem.  The  resulting  com¬ 
puter  program  was  then  used  to  evaluate  reflected  and  trans- 


I - 

mitted  waves  arising  from  a  spherical  wave  interacting  with 
a  finite-thickness  fluid  plane.  Calculations  were  performed 
for  a  variety  of  fluid  properties  and  for  test  frequencies  1  and 
5  kHz.  The  fluid  properties  chosen  do  not  correspond  to 
those  of  any  real  fluids,  but  rather  are  representative  of  the 
densities  and  sound  speeds  characteristic  of  the  solids  of 
interest  (i.e.,  metals  and  elastomers).  The  density  range  of 
the  fictitious  fluids  considered  was  0.5-0. 8  g/cm\  The 
sound-speed  range  considered  was  1.5  X  105-7.0X  105  cm/s. 

The  results  were  next  compared  to  numerical  solutions 
obtained  using  a  much  simpler  computer  program  that  im¬ 
plements  the  solution  to  the  associated  problem  involving 
the  interaction  of  a  normally  incident  plane  wave  with  a  fluid 
plane.  The  comparison  was  achieved  by  first  computing  the 
reflected-wave  pressure  at  the  point  r  =  0,  z  =  0  for  the  case 
of  the  incident  spherical  wave.  The  transmitted  wave  pres¬ 
sure  at  the  point  r  =  0,  z  =  —  /  was  also  computed  for  this 
case.  The  results  so  obtained  were  compared  to  the  corre¬ 
sponding  results  obtained  using  plane- wave  theory,  evaluat- 
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ed  at  the  relevant  surfaces.  In  all  cases  considered,  the 
spherical-wave  solution  differed  very  little  from  the  plane- 
wave  solution.  For  example,  consider  an  incident  spherical 
wave  of  1-kHz  frequency,  originating  from  a  source  located 
at  z'  =  200  cm,  interacting  with  the  fictitious  fluid  having 
density  1.2  g/cmJ.  sound  speed  3  X  10s  cm/s,  and  1-cm  thick¬ 
ness.  In  this  case,  the  spherical-wave  reflection  coefficient, 
evaluated  on  the  surface  of  the  fluid  plane  that  is  facing  the 
source  and  normalized  to  the  incident  wave  at  the  origin  of 
coordinates,  is  9.33  X  10“4-2.35x  10“2y,  while  that  for  the 
plane  wave  is  6. 12  X  10“4-2.08X  10 “7.  The  spherical-wave 
transmission  coefficient  for  this  case,  evaluated  on  the  sur¬ 
face  of  the  fluid  plane  that  is  opposite  the  source  and  also 
normalized  to  the  incident  wave  at  the  origin  of  coordinates, 
is.  0.993  +  2.97x10“ 7,  while  that  for  the  plane  wave  is 
0.994  +  2.93  X  10“  7. 

Next,  the  computer  program  that  implements  the  solu¬ 
tion  to  the  full  elasticity  problem  was  used  to  evaluate  the 
reflected  and  transmitted  waves  arising  from  the  interaction 
of  a  spherical  wave  with  a  finite-thickness  fluid  plane.  In 
order  to  do  this,  an  extremely  small  value  ( 10“A  dyn/cm2) 
for  the  shear  modulus  fi  was  used.  In  all  cases  tested,  the 
results  obtained  using  the  general  computer  program  agreed 
with  those  obtained  using  the  specialized  computer  program 
for  this  problem  to  within  the  precision  of  implementation. 

Taking  the  above  consistencies  as  indicating  that  the 
computer  program  that  implements  the  general  solution 
does  not  contain  simple  programming  errors,  cases  were 
next  considered  which  involved  a  solid  plate.  The  results  of 
several  representative  calculations  are  given  in  Table  I.  The 
first  column  of  this  table  gives  the  material  considered.  The 
second  column  gives  the  frequency  considered.  The  third 
through  sixth  columns  present  reflection  coefficients  (R) 
and  transmission  coefficients  ( 7~)  for  both  the  spherical-wave 
and  plane-wave  cases.  (The  coefficients  are  normalized  to 
the  incident  wave  at  the  origin  of  coordinates. )  The  reflec¬ 
tion  coefficients  are  evaluated  at  r  =  0,  z  =  0  and  the  trans¬ 
mission  coefficients  are  evaluated  at  r  =  0,  z  =  —  1  cm. 
Plate  thickness  in  all  cases  considered  is  /  =  1  cm.  Spheri¬ 
cal-wave  source  offset  is  z'  =  200  cm  in  all  cases.  The  shear- 


wave  speed  for  each  material  is  listed  for  reference  in  the  last 
column  of  the  table.16 

Note  that  for  several  of  the  spherical-wave  cases  pre¬ 
sented,  the  amplitude  of  the  reflected  and/or  transmitted 
waves  exceeds  that  of  the  incident  wave  evaluated  at  the 
origin  of  coordinates.  Although  this  unexpected  result 
would  be  indicative  of  error  in  a  normal-incidence  plane- 
wave  calculation,  *his  is  not  necessarily  so  for  an  incident 
spherical  wave,  since  the  reflected  and  transmitted  pressure 
amplitudes  will  each  vary  with  r  in  this  case.  Nonetheless, 
one  is  reluctant  to  accept  the  validity  of  such  computed 
“overpressures”  without  further  investigation. 

In  order  to  determine  whether  the  unexpected  results  of 
Table  I  might  have  possibly  arisen  from  either  numerical 
errors  or  computer  programming  errors,  an  energy -conser¬ 
vation  analysis  was  performed.  This  calculation  involved  nu¬ 
merically  performing  the  power  integral17 

E  =  —  Re  f  p*(r)v(r)-«(2nr dr).  (31) 

2  Jo 

The  asterisk  in  Eq.  (31)  denotes  complex  conjugate.  This 
integral  was  evaluated  over  both  the  front  and  back  surfaces 
of  the  plate,  and  the  result  was  compared  to  the  total  power 
incident  on  the  infinite  plate.  (The  total  incident  power  is 
equal  to  one-half  the  total  power  radiated  by  the  point 
source. )  The  particle  velocity  v  required  in  Eq.  (31)  was  eval¬ 
uated  numerically  at  each  point  of  integration  by  using  the 
gradient  of  the  appropriate  pressure  term.  [In  essence,  the 
partial  derivative  with  respect  to  z  of  Eq.  (17)  was  used  to 
compute  the  required  component  of  the  particle  velocity  as¬ 
sociated  with  the  reflected  pressure,  and  this  same  deriva¬ 
tive  of  Eq.  (18)  was  used  to  compute  the  required  component 
of  the  particle  velocity  associated  with  the  transmitted  pres¬ 
sure.  The  exact  relationship  between  each  pressure  gradient 
and  the  relevant  particle  velocity  is  provided  by  Euler's 
equation.  ]  The  integral  of  Eq.  (31)  was  performed  by  using  a 
32-point  Gauss-Legendre  quadrature  scheme. 12  Each  of  the 
required  pressures  and  velocities  was  evaluated  at  each  of 
the  32  quadrature  points.  Since  each  required  velocity  and 
pressure  is  defined  in  terms  of  an  integral,  each  power  calcu- 


TABLE  l.  Spherical-wave  complex  reflection  coefficients  R  and  transmission  coefficients  7  compared  to  plane-wave  coefficients  for  a  variety  of  materials, 
each  of  /  =  I -cm  thickness.  PMM  denotes  polymethylmethacrylate.  The  fluid  medium  is  water.  Spherical-wave  source  offset  z  in  each  case  is  200  cm  R  is 
evaluated  at  the  origin  of  coordinates.  T  is  evaluated  at  the  point  (0,0,  —  I) .  Each  coefficient  is  normalized  to  the  incident  w  ave  evaluated  at  the  origin  of 
coordinates.  ( See  Fig.  1 . ) 


Frequency 

(kHz) 

R 

r 

Shear  speed 

(cm/s) 

Spherical  wave 

Plane  wave 

Spherical  cave 

Plane  wave 

Aluminum 

-  0.482  +  0.140; 

0.003  -  0.056; 

1.50-0.223; 

0.997  +  0.057; 

2.98  X  10' 

Brass 

-  0.078  +  0.204; 

0.031  -0.172; 

1.02  -  0.350; 

0.969  +  0.173; 

2.11  X  10' 

Lead 

1 

0.104-0.185; 

0.053  -  0.223 j 

0.896  +  0.185; 

0.947  +  0.225; 

6.98  >.  104 

PMM 

-0.174  +  1.15; 

0.001  -  0.019; 

0.792  +1.18; 

0.999  +  0.030; 

1.34x10' 

Steel 

-0.616-1.07; 

0.025-0.157; 

1.69  F  1.02; 

0.974  +  0.157; 

3.28x  10' 

Aluminum 

-0.133-0.063; 

0.074  -  0.260; 

0.355  +  0.562; 

0.925  +  0.264; 

2.98  x  10' 

Brass 

1.17-0.614; 

0.442  -  0,495; 

1.22-0.002; 

0.558  +  0.498; 

2.11x10' 

Lead 

5 

0.575  -0,457 

0.584  -  0.488; 

0.420  +  0.468; 

0.415  +  0.498; 

6.98  xlO4 

PMM 

0.096  -  0.237; 

0.014  -  0.094; 

1.03  -  0.001; 

0.983  +  0.149; 

1.34  x  10' 

Steel 

0.625  +  0.936; 

0.394  -  0.049; 

0.015  -  0.844; 

0.605  +  0.489; 

3.28  x  10' 
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lation  thus  involved  6S  numerical  integrations  ( 32  velocity 
integrals,  32  pressure  integrals,  and  1  integral  for  the  calcu¬ 
lation  of  the  power).  Of  course,  the  integral  of  Eq.  (31) 
should  actually  be  performed  for  the  upper  limit  r'  —  oo . 
However,  if  a  sufficiently  large  value  for  r  is  chosen,  the 
vast  majority  of  the  power  will  be  included. 

Table  II  presents  results  of  the  power  calculations  for 
the  1-kHz  case  involving  a  steel  plate  of  1-cm  thickness. 
Note  from  the  spherical-wave  values  listed  in  Table  I  that 
the  pressure  amplitude  of  the  reflected  wave  exceeds  that  of 
the  incident  wave  by  more  than  23%  in  this  case.  The  trans- 
mitted-wave  pressure  amplitude  is  nearly  twice  that  of  the 
incident  wave.  Hence,  this  is  a  particularly  severe  case  to 
examine. 

In  Table  II,  the  first  column  gives  the  upper  limit  of 
integration  r'  used  in  evaluating  the  integral  of  Eq.  (31).  It 
represents  the  radius  of  a  circle  drawn  on  either  side  of  the 
plate  having  its  center  either  (i)  at  the  origin  of  coordinates 
( on  the  side  of  the  plate  facing  the  source )  or  ( ii )  at  the  point 
(0,0,  —  /)  (on  the  side  of  the  plate  opposite  the  source).  One 
such  circle  may  be  taken  to  be  that  illustrated  in  Fig.  1,  if  the 
r  coordinate  is  taken  to  be  equal  to  r' .  For  the  purposes  of  the 
following  discussion,  we  will  refer  to  either  of  these  circles 
as  a  “circle  of  observation." 

The  quantity  E,  (r  )  represents  the  power  contained  in 
that  portion  of  the  incident  wave  which  strikes  the  infinite 
plate  within  the  circle  of  observation  depicted  in  Fig.  1.  The 
quantity  Er  (r)  denotes  the  power  contained  in  that  portion 
of  the  reflected  wave  that  also  passes  through  the  circle  of 
observation  on  the  side  of  the  plate  facing  the  source.  Simi¬ 
larly,  the  quantity  E,(r')  denotes  the  power  contained  in 
that  portion  of  the  transmitted  wave  that  passes  through  the 
circle  of  observation  on  the  side  of  the  plate  that  is  opposite 
the  source.  All  values  are  expressed  as  a  percentage  of  the 
total  power  incident  on  the  infinite  plate.  The  power  balance 
can  be  determined  by  summing  the  values  in  the  Er  (/•')  and 
E,(r )  columns,  and  comparing  the  result  to  100%  (the  total 
incident  power).  The  sum  Er(r')  -f  E,(r')  is  given  in  the  last 
column  of  Table  II.  As  can  be  seen,  the  total  never  exceeds 
100%.  The  10000-crr  radius  considered  in  Table  II  is  the 
greatest  value  that  could  be  accommodated  by  the  numerical 
scheme  used  to  evaluate  the  required  integrals.  It  should  be 
noted  that  this  radius  accounts  for  98.4%  of  the  total  inci¬ 
dent  power. 

By  comparing  the  values  in  the  E,  (r1)  column  of  the 

TABLE  II.  Power  balance  for  a  1-kHz  spherical  wave  originating  from  a 
source  located  at  z  =  200 cm  and  interacting  with  a  steel  plate  of  1  .cm  thick¬ 
ness  and  infinite  lateral  extent.  £, (/) — power  incident  on  the  relevant 
“circle  of  observation"  ( radius  r’);  E,  ( / )  — power  reflected  through  the  rel¬ 
evant  circle  of  observation;  E,  (ri) — power  transmitted  through  the  rel¬ 
evant  circle  of  observation.  Results  are  expressed  as  a  percentage  of  the 
total  power  incident  on  the  infinite  plate. 


r’(cm) 

£,(/) 

EAn 

£,(/■') 

E,{r)  +£,(/) 

100 

10.5 

4.18 

17.1 

21.3 

250 

37.5 

15.8 

24.2 

40.0 

1000 

80.3 

51.4 

31.0 

82.4 

10000 

98.4 

62.8 

37.0 

99.8 
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table  with  the  values  in  the  column  that  gives  the  sum  of 
E,  ( r‘ ’ )  and  E,  ( r' ) ,  it  becomes  clear  that  the  power  required 
to  create  the  computed  overpressures  is  provided  by  power 
contained  in  portions  of  the  incident  spherical  wave  that 
strike  the  plate  in  regions  that  are  outside  the  circle  of  obser¬ 
vation.  This  conclusion  follows  from  the  fact  that  the  sum 
E,  (/)  +  E,  (r‘)  is  always  greater  than  E,  (r‘),  although  the 
difference  {[Er(r')  -f  E,  (r‘)  ]  —  E,(r')}  clearly  appears  to 
approach  zero  as  r'~.  oc.  Thus  the  additional  power  arrives 
at  the  observation  point  either  by  propagating,  within  the 
plate  material,  from  regions  outside  of  this  circle  to  regions 
interior  to  it,  and  then  being  reradiated  into  the  fluid,  or  by 
direct  radiation  into  the  fluid  by  regions  of  the  plate  that  are 
outside  of  this  circle.  Of  course,  both  of  these  alternatives 
may  contribute  to  the  phenomenon.  These  unexpected  re¬ 
sults  appear  to  be  a  consequence  of  solving  the  full  three- 
dimensional  elasticity  theory.  Several  attempts  which  were 
considered  to  explain  these  overpressures  using  simplified 
theories  were  unsuccessful. IK  W  (For  example,20  the  flex¬ 
ural-wave  coincidence  frequency  for  the  steel  case  consid¬ 
ered  is  approximately  23.7  kHz. ) 

We  can  gain  further  understanding  of  how  it  is  possible 
for  overpressures  to  be  present  and  yet  for  energy  to  none¬ 
theless  be  conserved  by  investigating  the  radial  dependence 
of  the  reflected  and  transmitted  pressures.  Three  cases  are 
considered  in  Figs.  2  and  3. 

Figure  2  ( a )  -  ( c )  presents  reflected-wave  pressure  ampli- 
tudes  as  a  function  of  r,  evaluated  at  z  =  0  and  z  —  150  cm, 
for  steel,  polymethylmethacrylate,  and  lead  samples  of  1  -cm 
thickness,  respectively.  Soiree  offset  is  z'  =  200  cm  and 
source  frequency  is  1  kHz.  Figure  3  ( a )  -  ( c )  presents  transmit  - 
ted-wave  results  for  these  same  cases,  except  that  the  pat¬ 
terns  are  evaluated  at  z  =  —  lcmandz=  —  151cm. Ascan 
be  seen,  each  of  these  cases  has  an  associated  beam  pattern. 
Therefore,  although  the  reflected  and/or  transmitted  wave 
pressure  amplitudes  can  exceed  that  of  the  incident  wave  in 
the  vicinity  of  r  =  0,  these  pressure  amplitudes  eventually 
drop  below  that  of  the  incident  wave  at  greater  values  of  r. 
Thus  it  is  not  unreasonable  that  the  integrated  power  over 
both  sides  of  the  plate  can  be  equal  to  the  total  incident 
pow  r,  even  though  the  reflected  and/or  transmitted  wave 
pressure  amplitudes  exceed  that  of  the  incident  wave  in  cer- 
ta;n  limited  regions. 

Taking  this  power  calculation  as  sufficient  for  establish¬ 
ing  the  validity  of  the  present  results,  we  next  investigate  a 
possible  cause  of  the  computed  overpressures.  The  results 
presented  in  Table  I  tend  to  support  the  notion  that  the  effect 
is  associated  with  the  nonzero  shear-wave  speed  associated 
with  solids.  (A  similar  conclusion  was  reached  in  Ref.  1 
concerning  a  cylindrical  scatterer. )  Evidence  for  this  is  pro¬ 
vided  by  first  noting  the  small  difference  between  the  spheri¬ 
cal-wave  and  plane-wave  results  presented  for  lead  in  Table 
I.  (Note  that  the  shear  speed  of  lead  is  about  one-half  the 
sound  speed  in  water. )  Since  lead  has  a  relatively  low  shear 
speed,  this  minor  difference  between  the  spherical-wave  re¬ 
sults  and  the  plane-wave  results  may  be  primarily  attribut¬ 
able  to  the  simple  phase  and  amplitude  variation  of  the  inci¬ 
dent  wave  across  the  plate  surface,  rather  than  being  due  to  a 
significant  shear  effect.  (Recall  also  the  close  agreement 
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FIG.  2.  Reflected-wave  pressure  amplitude,  as  a  function  of  r,  resulting 
from  a  1-kHz  incident  spherical  wave  originating  from  a  source  located  at 
z'  —  200  cm.  Solid  line:  z  =  0;  dashed  line:  z  —  150  cm;  sample  thickness:  1 
cm.  Incident  wave  at  r  =  O.z  =  0  is  taken  to  be  of  amplitude  1.0:  (a)  steel;  (b) 
polymethylmethacrylate;  (c)  lead. 


FIG.  3.  Transmitted-wave  pressure  amplitude,  as  a  function  of  r,  resulting 
from  a  1-kHz  incident  spherical  wave  originating  from  a  source  located  at 
z  =  200cm.  Solid  line:z  —  —  lcm;dashedline:z  =  —  151cm;samplethick- 
ness:  1  cm.  Incident  wave  at  r  —  0,  z  =  Ois  taken  to  be  of  amplitude  1.0:  (a) 
steel;  (b)  polymethylmethacrylate;  (c)  lead. 
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between  spherical-wave  theory  and  plane-wave  theory  for 
the  case  of  the  fluid  scatterer  discussed  earlier. )  Note,  how¬ 
ever,  for  the  materials  having  shear  speed  comparable  to  or 
significantly  greater  than  the  sound  speed  in  water,  the  dif¬ 


ferences  between  plane-wave  and  spherical-wave  results  are 
much  more  pronounced. 

We  further  investigate  the  notion  that  the  computed 


1291 


J  Acoust.  Soc.  Am.,  Vol.  83,  No.  4,  April  1988 


Jean  C.  Piquette:  Spherical-wave  scattering  by  a  plate 


1291 


overpressures  are  associated  with  the  influence  of  shear  ef¬ 
fects  in  Fig.  4(a)  and  (b).  In  these  figures,  the  reflected  and 
transmitted  wave  pressure  amplitudes,  evaluated  on  the  rel¬ 
evant  surfaces  of  the  plate,  are  evaluated  for  a  continuum  of 
fictitious  materials.  As  with  the  previous  figures,  the  pres¬ 
sure  amplitudes  are  normalized  to  that  of  that  incident  wave 
at  the  origin  of  coordinates.  Each  of  these  fictitious  materi¬ 
als  is  taken  to  have  the  same  density  and  the  same  longitudi¬ 
nal  wave  speed  as  that  of  steel.  However,  the  shear  modulus 
is  allowed  to  vary  between  zero,  at  the  left  end  of  the  abscis¬ 
sa,  and  the  true  shear  modulus  of  steel,  at  the  right  end  of  the 
abscissa.  (The  shear  modulus  is  normalized  to  that  of  steel, 
so  that  the  abscissa  runs  from  0  to  1.)  Note  that  the 
pressure  amplitude  corresponding  to  u/ps  =  0  in  each  of 
these  figures  is  less  than  the  incident  wave  pressure  ampli¬ 
tude  at  the  origin  of  coordinates,  but  that  this  no  longer 
holds  true  as p/fi,  departs  from  zero.  It  is  interesting  to  note 


FIG.  4.  Reflected-  and  transmitted-wave  pressure  amplitudes  arising  from 
an  incident  spherical  wave,  as  a  function  of  normalized  shear  modulus,  for  a 
continuum  of  fictitious  quasisteel  samples.  The  quantity  fit  denotes  the 
shear  modulus  of  steel.  For  all  cases,  the  longitudinal  wave  speed  in  the 
sample  is  held  fixed  at  the  true  value  for  steel.  Thus ////<,  =  0.0  corresponds 
to  a  fictitious  fluid  sample  having  the  same  sound  speed  as  the  longitudinal 
wave  speed  in  actual  steel,  andfi/fi,  =  1 .0  corresponds  to  actual  steel.  Sam¬ 
ple  thickness:  1  cm;  source  offset:  200  cm;  frequency:  1  kHz:  (a)  reflected- 
wave  pressure  amplitude,  evaluated  at  r  =  0,  z  =  0;  (b)  transmitted-wave 
pressure  amplitude,  evaluated  at  r  =  0,  z  =  —  1  cm. 
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that  the  departure  of  the  results  from  those  at  p/p,  =  0  is 
nonmonotonic,  and  that  there  are  significant  intervals  for 
which  the  computed  pressures  are  significantly  less  than 
those  expected  from  plane-wave  theory.  ( The  values  present¬ 
ed  for  p/p,  =  0  do  not  differ  significantly  from  the  plane- 
wave  results. ) 

III.  SOME  CONSEQUENCES  FOR  PANEL 
MEASUREMENTS 

When  an  unexpected  result  such  as  that  reported  here  is 
obtained,  one  is  always  concerned  that  the  “phenomenon” 
may  be  a  consequence  of  mathematical  errors  or  computer 
programming  errors.  However,  it  would  be  rather  surprising 
if  such  errors  yielded  results  that  are  consistent  with  the  law 
of  conservation  of  energy.  Also,  the  same  computer  program 
that  predicts  the  unexpected  results  also  predicts  very  rea¬ 
sonable  results  for  the  cases  in  which  the  shear  speed  in  the 
solid  is  low  (e.g.,  lead).  Nevertheless,  these  results  are  suffi¬ 
ciently  surprising  that,  even  with  the  validity  test  provided 
by  the  energy  conservation  law,  one  still  tends  to  be  skepti¬ 
cal. 

An  attempt  was  made  to  acquire  experimental  data  to 
verify  the  theory.  Specifically,  reflection  and  transmission 
data  were  acquired  for  the  frequency  range  1-5  kHz  using 
both  steel  and  PMM  samples.  The  samples  were  each 
square-shaped,  76  cm  on  a  side.  The  experimental  param¬ 
eters  were  the  same  as  those  used  to  obtain  the  theoretical 
results  presented  earlier. 

These  measurements  were  found  to  be  in  reasonably 
good  agreement  with  simple  plane-wave  theory  and.  hence, 
in  apparent  disagreement  with  the  predictions  of  Table  I. 
Insight  into  the  explanation  for  these  negative  results  can  be 
gained  by  again  referring  to  Table  II.  Note  that,  for  the 
r  =  100-cm  radius  case,  the  total  reflected  -I-  transmitted 
power  is  21.3%  of  the  total  incident  power.  However,  the 
E,  ( r' )  column  shows  that  only  10.5%  of  the  total  incident 
power  actually  strikes  the  relevant  circle  of  observation  for 
this  case.  This  means  that  if  a  measurement  were  carried  out 
using  a  circularly  shaped  steel  plate  of  100-cm  radius  and  1- 
cm  thickness,  quantitative  agreement  with  the  predicted  re¬ 
sults  could  not  possibly  occur  since,  as  previously  men¬ 
tioned,  much  of  the  power  required  to  create  the  predicted 
overpressures  originates  from  portions  of  the  incident  wave 
that  would  strike  the  associated  infinite  plate  outside  the 
circle  of  observation.  For  quantitative  agreement  with  the 
theory,  a  test  plate  of  sufficiently  large  radius  would  be  re¬ 
quired  that  the  total  power  actually  incident  on  the  test  plate 
would  not  differ  significantly  from  the  total  reflated  plus 
transmitted  power  passing  through  the  associated  circles  of 
observation.  [The  250-cm  radius  case  considered  in  Table  II 
may  be  a  good  candidate,  since  the  sum  Er(r')  +  E,(r')  ex¬ 
ceeds/!,  (r")  by  only  2.5%  of  the  total  incident  power  in  this 
instance.  ]  Thus  the  negative  results  observed  in  the  measure¬ 
ments  described  here  (which  were  made  using  a  square¬ 
shaped  sample  having  an  area21  less  than  20%  of  that  of  the 
100-cm  observation  circle  of  Table  II )  are  most  likely  a  con¬ 
sequence  of  the  very  small  size  of  the  test  sample,  and  of  the 
data  acquisition  technique.22  However,  these  negative  ex¬ 
perimental  results  and  the  theoretical  predictions  of  Table  II 

Jean  C.  Piquette:  Sphencal-wave  scattering  by  a  plate 


1292 


1292 


combine  to  yield  an  important  conclusion  regarding  panel 
test  measurements:  Measurements  made  on  test  panels  of 
large  cross  section  may  be  significantly  corrupted  by  wave- 
front  curvature  effects. 

The  larger  the  test  sample  used,  the  greater  the  contri¬ 
butions  to  the  observed  pressures  from  panel  regions  that  are 
far  from  the  receiver.  It  is  clear  from  Table  II  that  the 
greater  the  contributions  become  from  distant  panel  ele¬ 
ments,  the  greater  is  the  possibility  that  undesirable  wave- 
front  curvature  effects  will  corrupt  the  results.  It  is  general¬ 
ly  desired  to  obtain  the  plane-wave  response  of  a  test  panel. 
It  is  also  generally  believed  that  panels  of  large  cross  section 
are  superior  to  panels  of  small  cross  section  for  conducting 
panel  tests,  since  the  undesirable  influence  of  panel  edge 
diffraction  can  be  reduced  by  increasing  the  panel’s  lateral 
dimensions.  However,  it  is  clear  from  the  results  presented 
here  that  the  edge-diffraction  reduction  benefits  derived 
from  increasing  sample  size  may  be  more  than  offset  by  the 
detrimental  influences  of  increased  wave-front  curvature  ef¬ 
fects.  Of  course,  this  conclusion  is  based  on  the  assumption 
that  the  usual  practice  of  establishing  plane-wave  conditions 
by  using  the  farfield  of  a  source  is  employed.  If  a  suitable 
array  technique  is  used  instead,  the  spherical-wave  effects 
described  here  may  be  suppressed. 

It  should  be  noted  that  the  large  overpressures  predicted 
here  have  been  only  theoretically  observed  to  occur  for  ma¬ 
terials  having  a  shear  speed  comparable  to  or  greater  than 
the  sound  speed  in  the  surrounding  fluid  and  that  panels  are 
generally  fabricated  from  soft,  rubbery  materials.  It  should 
also  be  noted  that,  even  for  stiff  materials,  if  a  sufficiently 
large  loss  is  present,  these  overpressures  disappear.2'  On  the 
other  hand,  test  panels  are  frequently  fabricated  using  a 
steel  backing  plate  for  support.  It  may  very  well  be  that,  for 
sufficiently  large  samples,  this  steel  backing  plate  intro¬ 
duces  undesirable  spherical-wave  effects  into  the  observa¬ 
tions.24 

IV.  SUMMARY  AND  CONCLUSIONS 

The  problem  of  the  interaction  of  a  spherical  wave  with 
a  homogeneous  and  isotropic  solid  plate  of  infinite  cross  sec¬ 
tion  but  finite  thickness  was  considered  theoretically.  Re¬ 
sults  of  numerical  calculations  of  reflected  and  transmitted 
pressure  waves  for  a  variety  of  test  materials  immersed  in  a 
water  medium  were  presented.  The  theory  predicts  that 
when  the  shear  speed  in  the  solid  is  comparable  to  or  exceeds 
the  sound  speed  in  the  surrounding  fluid,  the  reflected-  and 
transmitted-wave  pressure  amplitudes  can  exceed  the  maxi¬ 
mum  value  of  the  incident-wave  pressure  amplitude  on  the 
surface  of  the  plate.  These  results  are  of  practir  fl  impor¬ 
tance  to  those  concerned  with  acoustical  measurements 
made  on  planar  structures.  However,  a  complete  theoretical 
interpretation  of  the  causes  of  the  results  was  not  presented 
herein,  and  this  awaits  further  investigation. 

An  experiment  was  performed  on  polymethylmethacry¬ 
late  and  steel  samples  using  experimental  parameters  con¬ 
sidered  in  the  theoretical  calculations.  The  measurements 
were  found  to  be  in  better  agreement  with  simple  plane-wave 
theory  than  with  the  spherical-wave  calculations  presented 
here.  These  negative  results  am  very  likely  a  consequence  of 


the  very  small  cross-sectional  area  of  the  test  samples  used. 
An  important  conclusion  of  the  present  work  is  that  mea¬ 
surements  made  on  test  panels  of  large  cross-sectional  area 
may  be  significantly  corrupted  by  wave-front  curvature  ef¬ 
fects.  The  practice  of  using  a  steel  backing  plate  to  support 
test  panels  may  also  introduce  undesirable  spherical-wave 
effects  into  the  observations,  particularly  for  samples  of 
large  ernss  section. 
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area  element  lirrdr  is  appropriate  for  angularly  independent  acoustic 
fields.  The  unit  vector  h  is  equal  to  k  on  the  surface  of  the  plane  that  faces 
the  source,  and  is  equal  to  —  k  on  the  surface  of  the  plane  that  is  opposite 
the  source. 

'“The  mechanisms  considered  were  thickness  resonances  involving  longitu¬ 
dinal  or  shear  waves,  critical  angle  phenomena,  flexural  wave  coinci¬ 
dences,  and  “refraction  arrival"  phenomena  (see  Ref.  2.  pp.  195-198). 
'"Another  possible  mechanism,  which  was  not  initially  considered  but  was 
suggested  by  a  reviewer  of  this  article,  is  the  generation  of  the  symmetric 
Lamb  wave  by  the  incident  wave  It  is  possible  that  such  waves  may  in¬ 
clude  a  component  directed  toward  the  origin.  These  w  aves  would  be  fo¬ 
cused  at  the  origin,  and  may  therefore  serve  to  explain  the  reported 
overpressures.  However,  a  thorough  investigation  of  this  possibility  is  be¬ 
yond  the  scope  and  intent  of  the  present  article.  For  a  discussion  of  the 
physics  of  the  generation  of  both  symmetric  and  antisymmetric  Lamb 
waves  in  an  elastic  plate  under  plane-wave  excitation,  the  interested  read¬ 
er  is  directed  to  R.  Fiorito,  W  Madigosky,  and  H.  Uberall.  “Resonance 
theory  of  acoustic  waves  interacting  with  an  elastic  plate,"  J  Acoust.  Soc. 
Am.  66.  1857-1866  (1979). 

"See,  for  example.  Ref.  17,  p  128.  Eq.  (3-5.8),  and  the  discussion  in  the 
paragraph  following  this  equation.  Although  this  equation  is  strictly  valid 
only  for  thin  plates,  the  computed  coincidence  frequencies  are  sufficiently 
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different  from  the  frequencies  of  interest  that  this  consideration  is  not 
important  here. 

2lThe  tendency  exhibited  by  the  values  n'e**"*-d  in  Table  II  is  for  the  total 
power  propagating  away  from  the  plate  through  the  circles  of  observation, 
E,(r' )  +£,(/),  to  exceed  that  which  is  incident  on  the  circle  of  observa¬ 
tion  that  is  facing  the  source,  E,  (r  ),  by  an  increasing  percentage  for  de¬ 
creasing  values  of  r‘ .  Thus,  the  smaller  the  tc„t  sample  used,  the  smaller  is 
the  likelihood  that  the  predicted  overpressures  will  be  observed. 

:  The  data  acquisition  technique  involved  the  use  of  pulsed  interrogating 
waves  and  the  use  of  an  observation  “gate.”  This  gate  eliminated  contribu¬ 


tions  to  the  observed  waves  originating  from  the  sample  edges. 

The  results  presented  in  Table  1  for  polymethylmethacrylate  include  loss. 
Thus  the  loss  factor  of  this  material  is  insufficient  to  eliminate  the  predict¬ 
ed  phenomenon,  given  a  sufficiently  large  sample.  It  should  also  be  noted 
that,  even  though  the  overpressures  disappear  when  a  sufficiently  large 
loss  is  iutioduced,  the  significant  differences  between  plane-wave  and 
spherical-wave  theory  persist. 

2<Of  course,  the  damping  influence  of  the  lossy  materials  affixed  to  the  steel 
backing  plate  may  inhibit  the  spherical-wave  effects.  This  is  currently  an 
unresolved  issue  that  must  be  answered  by  further  research. 
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cTrfcqlar  segment  [Fig  !%d)  j,  and  this  effect  increases  with 
an  inc><ase  in  the  value  of^/r,,,  which  makes  the  segment 
flatter,  approaching  a  rectanglbsas  shown  in  Fig.  3  for  the 
case  of  three 'fkgimuthal  partitions  (sot  —  3). 

iii.  concluding  Remarks 

The  foregoing  study  indicates  the  follow'ihg' 

( 1 )  One  azimuthal  partitibn.  does  not  raiseTty?  cutoff 
frequency  (the  smallest  cut-on  frequency)  at  all. 

(2)  Azimuthal  partitioning  can  at  best  (for  thp»e^rair 

muthal  partitions)  raise  the  cutoff  frequenjjKfJarame^r' 
{krmn  r,i)  from  1.84  to  3.83.  , 

(3)  Provision  of  a  radial  partiyatltm  addition,  >s  coun¬ 
terproductive  because  it  low-efstifu;  cutoff  frequency  param¬ 
eter  in  the  outer  annulapsCgments. 

The  last  observation  is  indeed  fortuitous  inasmuch  as  it 
would  be  extt^tifely  difficult  to  locate  loudspeaker.  in  the 
inner  segm*ms  of  the  duct  in  an  active  noise  control  system. ' 

Fipjflly,  it  may  be  pointed  out  that  in  the  foregoing  anal- 
ysis/fhe  wall  and  the  partitions,  Of  the  duct  have  been  as¬ 
signed  to  be  sufficiently  massive  and  rigid,  otherwise  there 
/would  be  sound-structure  interaction  altering  the  results  of 
the  investigation.  Besides  the  role  of  damping  has  also  been 


neglected  altogether.  In  any  gase  viseOus  damping  would 
have  only  a  marginal  effect  on Tji«'£ut-on  frequencies  of  high¬ 
er  order  modes. 
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Equation  (17)  was  erroneously  printed  twice  in  this  article,  while  Eq.  (18)  was  omitted.  Equation  (18)  should 
have  been 

„  _  r  b(0)J„{Pr)exp[jy  k  -  -  f37  {z' -z)] 
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